First observation of Kl — > 7T^e^ve~^e 
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This letter is the first report of the Kl — > tt^c^ ve'^ e~ decay. Based on 19208 ± 144 events, we 
determine the branching fraction, B{Kl n'^e^ve'^ e-;M,+ ,_ > 5 MeV/c\E*+^^ > 30 MeV) = 
(1.285±0.041)x 10"/ a,ndr{Ke3ee;M^+^- > 5MeY /c'^)/r{Ke3) = [4.57±0.04(stat) ±0.14(sysf)] x 
lO"'^. This ratio agrees with a theoretical prediction based on chiral perturbation theory (ChPT) 
calculated to (^(p*). The measured kinematical distributions agree with those predicted by just 
ChPT C(j)*), but show significant disagreement with ones predicted by leading order ChPT. 

PACS numbers: 13.20.-v, 13.20.Eb, 13.25.Es, 14.40.Aq, 12.39.Fe 



The semileptonic mode, Kl — > Tr^e^i^ (-^es) and 
its radiative mode, K^s-y, has been extensively studied 
In this letter, we introduce the semileptonic kaon 
decay mode Kl — > 7r^e^z/e+e~ (ifeSee)- We present the 
first measurement of its branching fraction, and the ratio 
of its decay width to that of the K^^ decay. As shown 
in Fig. [1] the Kl ^ ■K^e^ve^ e~ decay is dominated by 
Kes with a virtual photon, 7*, that converts internally 
into a real e+e~ pair. The amplitude of K^^-^ consists of 
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FIG. 1: Processes contributing to the Kl — > ■K^e^ve'^e~ de- 
cay: the virtual photon comes (a) from the charged pion, (b) 
from the K^z electron, and (c) from the kaon decay vertex. 
The process (c) includes the structure dependent (SD) ampli- 
tudes. The exchange diagrams have been omitted for clarity. 

two parts. One part is inner bremsstrahlung (IB) from 



the pion or the electron, which is illustrated by Fig.[T]-(a) 
and (b). The other part is the photon radiated from an 
intermediate hadronic state of the K-ir current, namely 
the structure dependent (SD) amplitude (ie. direct emis- 
sion) [ll,[3,|i], as illustrated by Fig.[T]-(c). 

The model describing the K-t: current is important for 
both studying K^z decays themselves and understanding 
low energy QCD. A powerful way to express the K-t: 
current is via Chiral Perturbation Theory (ChPT) 
ChPT has been developed based on the chiral symme- 
try of QCD, and it can be applied to all K^^ modes, 
including K^zee- In this letter, we compare our measure- 
ments against ChPT calculated to next-to-leading order, 
expanded to the fourth power of the momentum of chiral 
field p [NLO(/)]. 

We search for KeSee decays in data collected by the 
KTeV E799-II fixed target experiment, which ran at the 
Fermi National Accelerator Laboratory. Two parallel 
Kl beams were produced by 800 GeV/c protons from 
the Tevatron striking a BeO target. Following the tar- 
get were beam line elements to sweep away charged par- 
ticles, to absorb photons, and to allow for short-lived 
hadrons to decay away. The region from 95 m to 159 
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m downstream of the target was in vacuum, and defines 
the fiducial volume for Kl decays. Following a thin vac- 
uum window at the end of the decay region was a drift 
chamber spectrometer. The spectrometer had two pairs 
of drift chambers separated by an analysis magnet pro- 
viding a transverse momentum kick of 0.2 GeV/c. The 
momentum resolution of the spectrometer was measured 
to be (Jp/p — 0.016% X 0.38%, where p is the momen- 
tum of a charged particle in GeV/c. A set of transition 
radiation detectors (TRD) downstream of the spectrom- 
eter was used to distinguish pions and electrons. Farther 
downstream, there was a pure cesium iodide (Csl) elec- 
tromagnetic calorimeter, where the energy resolution for 
photons and electrons was (T£;/ii^ = 2%/\/E©0.45%, with 
E in GeV. Immediately upstream of the Csl calorimeter 
were scintillator hodoscopes, which served as the charged 
particle trigger. Behind the Csl was a 24 interaction- 
length steel filter and a set of scintillator hodoscopes 
to identify muons. Lead-scintillator counters were posi- 
tioned around the vacuum decay region, the spectrometer 
and the calorimeter, to reject events with particles escap- 
ing these detectors. We analyzed data acquired during 
the 1997 run. A detailed description for this experiment 
and analysis can be found in Ref . 0, 0] ■ 

The event reconstruction begins with the identification 
of four charged tracks coming from a vertex in the decay 
region. The charged tracks are identified as 7r^e^e+e~ 
using E /p, the energy reconstructed in the Csl calorime- 
ter divided by the momentum measured in the spectrom- 
eter. Pions tracks are required to have E/p less than 0.9, 
which selects 99.2% of all pions. Electron tracks are re- 
quired to have E/p between 0.93 and 1.15, and be tagged 
by the TRD system. The E/p and TRD requirements se- 
lect 95.0% of all electrons, while rejecting 99.95% of all 
pions. Since the ifeSee decay has three electrons, there 
are two candidates for the e+e~ pair. Although each 
event must include both amplitudes in which one of e^e~ 
pairs is from the virtual photon, in this letter, we define 
the pair which has smaller invariant mass as the "e^e^ 
pair" , and call the remaining electron as "e^^g" . The 
amplitude in which the smaller invariant mass e^e~ pair 
emerges from the virtual photon is greater in contribu- 
tion than the other amplitudes. Because the neutrino is 
not observed, there is a two-fold ambiguity for the par- 
ent kaon energy. The higher kaon energy solution was 
required to be less than 200 GeV. 

Monte Carlo (MC) simulations are used to understand 
the acceptance of the signal mode, background modes, 
and the normalization mode Kl tt+tt^tt^ where tt^ 
denotes the 7r° — > e"'"e~7 decay. To simulate Ke^ee, the 
matrix element was computed using ChPT[NLO(p'')] as 
described by Tsuji et alj^. Bremsstrahlung photons 
from four charged particles in i^TeSee are simulated us- 
ing the PHOTOS program, which includes the interference 
terms [l3| • The number of background MC events are 
generated according to their branching fractions and Kl 



flux estimated using the normalization mode. 
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FIG. 2: The fc+_o distributions for data and MC after all 

analysis requirements except for o'. The open histogram 

is the sum of all MC simulated modes (MC Total). The ver- 
tical line and arrow show the accepted region for the signal 
candidates (/fc+-o< -0.002 GeVVc^)- 



The backgrounds are reduced with a combination of 
particle identification, kinematical cuts, and vetos as in- 
dicated in Table [D The major background for the i^eSee 
mode is due to Kl — > 7r"'"7r~7r^. MC studies show that 
42% of the background arising from the Kl tt+tt^tt^ 
decays are caused by one of the pions misidcntificd as an 
electron. The rest are due to an external photon con- 
version in the detector material and a missing pion and 
electron. Missing pions are due to track hits corrupted by 
hadronic interactions in the detector material, while the 
analysis magnet causes low momentum tracks to escape 
the detector. We suppress the Kl ^ tt+tt^tt? 
ground using the kinematical variable. 



back- 



(1) 

where Mk and M^o are the kaon and tt^ masses, respec- 
tively. M^efces the invariant mass of the tt^ and ej^g, 
while the charged pion mass is assigned to the ej^g. pt 
is the transverse momentum of the Tr^eJ^g system. For 

Kl 7r+7r~7r° decays, fc_| o is the squared longitudinal 

momentum of the tt^ in the frame in which the momen- 
tum of TT+TT" system is transverse to the Kl direction, so 

that will be positive definite (Fig. [2|). On the other 

hand, for ifeSee events, o tends to have an unphysical 

value (fc^ o< 0). Therefore, we require k^ o < —0.002 

GeV^/c^ for the signal events. The radiative Ke3 decay 
with an external photon conversion are rejected by re- 
quiring Mg+g- > 5 MeV/c^. We also considered and 

and 



simulated backgrounds due to Kl 
Kl — » 7r+7r^7r4g. These backgrounds are small mainly 
because of the small branching fractions. 

In addition to the backgrounds in Table [H we also con- 
sidered and simulated two coincident Kl n'^e'^i' de- 
cays, and the S — > A(^ ptt~)tt^ decay; both of these 
backgrounds are negligible. After all cuts, we are left 
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TABLE I: The estimated backgrounds (BG) for the ifeSee analysis and the dependence on analysis cuts. The third column 
show the BG/Signal reduction factor stemming from the combined particle ID, fiducial, and veto cuts. The fourth column 
shows the same effect for the kinematic cut. 



Background 
decay mode 


BG/Signal 
before cuts 


BG/Signal 
after all cuts 


BG/Signal Reduction Factor 
Particle ID+Fiducial+Veto Cuts Kinematic cut 




117 


0.0182 


403 16 




0.048 


0.0167 


1.23 2.4 


Kl tt tt 7r4e 


0.307 


0.0096 


4.48 7.2 


Kl TT^e^wy, 7 conversion 


15.0" 


0.0081 


20.4 91 



"The minimum simulated photon energy was 1 keV. 



with a sample of 20225 events. The estimated total num- 
ber of background events after all cuts is 1017.1 ± 24.7, 
representing (5.03 ± 0.12)% of the signal sample. 

The normalization mode {K]^ tt+tt^tt^) were col- 
lected with the same conditions as the signal mode. They 

are analyzed using identical cuts, except that the 

requirement is reversed to be o> —0.002 GeY^/c^. 

We ignore the photon in the decay in order to make the 
analysis more similar to the signal mode analysis, which 
has a missing neutrino. The only significant background 



for the normalization analysis is the Kl 



de- 



cay followed by external conversion of one of the 7r° pho- 
tons. This background is determined by MC simulations 
to be (0.558 ± 0.005)% of 1250828 normahzation events 
after all cuts. 



The acceptance ratio of K^see to Kl 



depends on the efficiency ratios of an electron and a 
charged pion, since the signal mode has one more electron 
and one less pion compared to the normalization mode. 
Therefore, the efficiencies of an electron and a pion by 
the particle identification cuts are measured from data 
by tagging electrons and pions in Kl — > tt+tt^tt^ and 
Kl 7r"'"7r~7r°^ decays, respectively. The differences 
of the efficiencies between data and MC lead to a correc- 
tion applied to the MC signal acceptance, which becomes 
/e/, - 0.9955. 

After applying background subtractions and efficiency 
corrections, the decay width ratio of KeSee to Kl 
7r+7r~7r^ is 

7^(fce3ee/+ - Od) 

^ r(j^e3ee;M,+,- > 5MeV/c^g^^_ > 30MeV) 

r{KL tt+tt-ttOj) 
= [8.54 ± 0.07{stat) ± 0.13(sys<)] x 10"^ (2) 

where i?*+g- is the energy of e~^e~ pair system in the 
kaon rest frame. The acceptance of the signal events 
generated above the cut-off values is (0.8986±0.0025)%, 
and the acceptance of the normalization mode is 
(0.4947±0.0009)%. Table E] hsts the systematic errors 
in 7?.(fce3ee/H — O^i). The largest systematic error is the 
uncertainty in the number of the Kl 7r"'"7r~7r^ de- 
cays. The number of Kl — > tt+tt^tt^ decays measured 



TABLE II: Systematic uncertainties in the ratio of decay 
widths, TZ{ke3ee/-\ — Od), see Eq. [2]in the text. 



Source of 


Uncertainty on 


uncertainty 


7^(fce3ee/+ - Od)(%) 


Unobserved photon 




in normalization analysis 


±1.03 


Vertex cut 


±0.7 


Radiative corrections 


±0.51 


Corrections for the efficiency difference ±0.46 


Ek distribution 


±0.35 


Cut-off on the M^+^- 


-0.18 


Background estimations 


±0.05 


MC statistics 


±0.32 


Total systematic uncertainties 


±1.5 



using the photon (full reconstruction measurement) is 
(0.88±0.51)% smaller than the analysis ignoring the pho- 
ton. With this value and the systematic error in the 

we 



full reconstruction measurement of Kl — *■ tt 
assign a 1.03% systematic error on TZ{ke3ee/-\ — Od)- 
The second largest systematic error is based on a slight 
data-MC discrepancy in the distribution of the vertex x^, 
which indicates the quality of the four track vertex. The 
next largest systematic error is uncertainty in the treat- 
ment of radiative corrections. With inner bremsstrahlung 
photons generated in the MC by the PHOTOS program, 
the signal acceptance decreases by 3.6%. To confirm this 
acceptance loss, Kl — > ve'^ j (KeSeej) events are 
identified and compared between data and MC. We as- 
sign a systematic uncertainty from the error in the KeSecy 
measurement, although the number of KeSee-y evnents, 
935, is consistent with MC prediction. The probability to 
miss the tt track due to hadronic interactions in the TRD 



is determined by the GEANT program [ll| . The correction 
applied to the MC signal acceptance is 1.0328 ± 0.0045. 
We also estimate the uncertainties due to the E/p and 
TRD requirements. The total error in our estimate of 
the efficiency difference is ±0.46%. 

The KeSee branching fraction with statistical and sys- 
tematic uncertainty using B{Kl — > 7r+7r~7r°)= (12.56 ± 
0.05)% and B(7r° ^ 6+5-7)= (1.198 ± 0.032)% ^ is 

B{Ke3ee;M,+,- > 5MeV/c2,£;:+^_ > 30MeV) 
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[1.285 ± O.OlO(stat) ± OMO{syst)] x IO- 



CS) 



The systematic error is much larger than that of 
TZ{ke3ee/+ - Od), due to the 2.7% error on B{tt^ 

In the rest of this letter, we compare our results against 
the ChPT[NLO(p'')] description of the K-tt current. Us- 
ing the Ke3 branching fraction, B{Kes)= (40.53±0.15)% 
[121, we find 



n 



KeSee 



r(j^e3ee;M,+ ,- > 5MeV/c^) 

[4.57±0.04(stot)±0.14(syst)] x 10^^(4) 



The leading order ChPT and ChPT[NLO(p'^)] predic- 
tions for TZkeSee are (4.06 ± 0.11) x IQ-^, and (4.29 ± 
0.11) X 10~^ respectively. Although the experimental de- 
termination of TZKeSee iucludc all radiative effects, nei- 
ther the theoretical estimates of the numerator and de- 
nominator include these effects. To account for the lack 
of radiative treatment in the TZKeSee predictions, we as- 
sign an error of twice the value of SkcS- The variable 
^Kes is 0.013 ± 0.003, and parameterizes the increase in 
decay width of the K^^ mode due to radiative corrections 
[l3l|. The ChPT[NLO(p'^)] prediction is consistent with 
the measurement at the 1.6 cr level. 

As the K-TT form factor is parameterized by the square 
of the four momentum transfer to the leptons t = {pK — 
Ptt)^, higher order calculations of ChPT are sensitive to 
t. However, the K^^ee decay has a two-fold ambiguity in 
t due to the missing neutrino. To avoid this problem, we 
use the transverse momentum transfer as defined in [14| , 



Ml + Ml 



2M, 



K 



Ml 



(5) 



where M^r is the charged pion mass and px^^ is the trans- 
verse pion momentum. 

Figure [3] shows that the data ij^/M^ distribution 
agrees with the NLO(p'') calculation, but not with 
the leading order ChPT calculation. Figure [5] shows 
the invariant mass of the e+e~ pair, illustrating that 
ChPT[NL0(p4)] models well the ATeSee dynamics. The 
MT^eee,Meee, and Al-^e distributions are also well mod- 
eled with the ChPT[NLO(p'')] prediction. 

In summary, we find good agreement between our mea- 
surements and the NLO(p'*)ChPT calculation, while the 
leading-order ChPT calculation is disfavored. Finally, we 
note that Figure [3] is expected to receive contributions 
from both IB and SD amplitudes, with the IB ampli- 
tudes being dominant. The separation between IB and 
SD amplitudes has not been performed in the context of 
ChPT. Additional theoretical work is needed to extract 
the SD contribution. 

We gratefully acknowledge the support and effort of 
the Fermilab staff and the technical staffs of the partici- 
pating institutions for their vital contributions. We also 



1200 
-■1000 



ft 600 
•Z 400 



+ rdata , 




( 


a) 














h :MC 




800 


h :M< 
















ft- 600 














61.1/31 






1 -iff 
H 200 




43.4/31 


























1 2 


4 




7 "O I 




4 







oi (J+i+f*±t*-^ 



FIG. 3: Comparisons of the t±/Ml distributions for data 
(dots) and MC (histogram), (a) with MC-LO and (b) with 
MC-NLO(p''). The data-to-MC ratios at the bottom are fit 
to a straight line. 
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FIG. 4: Comparison of the M^+g- distribution for data 
(dots), and MC (histogram) with NLO(p*) correction. 
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